pss header will be provided by the publisher 



Comparison between models of insulator and semiconductor 
thin films islanding 

Francois Lallet* 1 , Alain Dauger** 1 et Nathalie Olivi-Tran*** 1 

1 SPCTS, UMR-CNRS 6638 47 a 73 Avenue Albert Thomas, 87065 Limoges cedex, FRANCE 



Mots-cle PVD-CVD, Sol-gel processing, Monte Carlo simulation, thin film islanding. 
PACS 04A25 

The synthesis of self-organized quantum dots (QD's) can be achieved through bottom up layer by layer 
deposition processes as chemical vapor deposition (CVD) or physical vapor deposition (PVD). However, 
QD's may also be synthesized via sol-gel route, which involves a spontaneous evolution from thin films to 
discrete QD's without further deposition. The aim of the paper is to discuss and compare the physical phe- 
nomena involved in QD's formation which initiate from thin film surface roughening between PVD-CVD 
and sol-gel synthesis approaches. We propose two simple physical models which are relevant to explain the 
fundamental differences between those methods. 
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1 Introduction 

Many authors have presented and studied epitaxial QD's growth through PVD-CVD synthesis pro- 
cesses. The physical phenomena associated with such approaches have been widely studied both theoreti- 
cally Q ED and experimentally H , in particular through the model system Ge/Si Q [10] flj] fl2) , 
because of its promising technological applications. However the synthesis of QD's epitaxially grown on 
a crystalline substrate can be achieved through a sol-gel approach. Indeed, Bachelet et al [13] have re- 
cently synthesized and studied the microstructure of zirconia QD's grown on a c-cut sapphire substrate 
during thermal annealing of a zirconia precursor xerogel thin film deposited on the substrate by sol-gel 
dip-coating. 

The aim of this article is to discuss the physical phenomena involved in these processes through the ana- 
lysis of two physical models, based on energetic considerations, which have been developed for numerical 
simulations. We suppose, for the sake of simplicity, that QD's synthesis is achieved without nucleation. 



2 Physical models 

We have developed a Monte Carlo (MC) algorithm to simulate the islanding of a thin film on a perfect 
crystalline subtrate during thermal annealing without further deposition [14|. The numerical thin film is 
divided into mesoscopic domains which are characterized by their height and their crystallographic orien- 
tation with regard to their nearest neighbors (NN) and the substrate. At each MC step, a domain i is chosen 
at random and the probability Pi to change its height (hi) and/or its crystallographic orientation (ci, di) 
is calculated through the classical Metropolis scheme [15 |. The energy of the domain i with regard to its 
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nearest neighbors is expressed as : 

, 2 \ NN , „ 2 \ NN j— NN 

Ei = li $>-c,)+72 )Y,^-d j )+Y{l+ V )J^—eY J (h l -h ] ) (1) 

\ I J i=1 J=l v 5 J=l 

where the first and second term of the right hand side of the equality correspond to the interfacial energy of 
the domain with regard to the substrate and to its neighbors respectively. The third term correspond to the 
surface energy related to the heights of the domains. 71 is the boundary surface tension (domain/domain), 
72 the interfacial surface tension (domain/substrate), I the distance between domain i and its nearest neigh- 
bors, Y the Young modulus, v the Poisson's ratio, D s the surface diffusion coefficient, 7 S the free surface 
tension, ks the Boltzmann constant and T the absolute temperature. 

We present another physical model which allows one to discuss, at least qualitatively, the main para- 
meters responsible for the morphological evolution of a thin film epitaxially grown on a perfect crystalline 
substrate through a deposition process (PVD or CVD). This model is inspired from the previous works of 
Kawamura [ 16] and Russo ifTTl who established MC algorithms at the atomic scale to model QD's forma- 
tion modes during deposition processes. The energy of an adatom i is computed as the sum of its bonding 
energy and elastic energy lfl6lllfT7l : 

E-i = NEg ~ {E w ith adatom i — ^without adatom i) (2) 

where N is the number of chemical bonds of the adatom i with its neighbors, Eb is the energy of a 
chemical bond and E the total elastic energy. From equation (2), we deduce a model in which an epitaxial 
monocrystalline thin film is deposited on a perfect crystalline substrate. We model this film as a cubic array 
of mesoscopic domains i of height hi and width i. Each domain is submitted both to its surface tension 
7 S and to an elastic stress field induced by the lattice mismatch between the film and the substrate, e. The 
volume energy density of a thin film of initial thickness h, due to free surface energy, is expressed as j s /h. 
Thus, the surface tension of a domain i can be written as (j s hi)/h and therefore the surface energy of 
domain i (related to NEb) induced by the free surface energy is {l 2 "f s hi)/h where £ 2 is the free surface 
area at the top of each domain. The elastic stress field leads to an elastic energy inside each domain of the 
film. This elastic energy (related to -E wml adatom i — -^without adatom i) can be expressed through linear elasticity 
theory. We suppose for the sake of simplicity that the elastic stress tensor is diagonal (for example for 
cubic phased materials). With the same assumption, e, which is the lattice mismatch between the film and 
the substrate, is constant in the horizontal plane. Thus the resulting force in the horizontal plane reduces 
to 2Ye£hi with £hi the area of surfaces of domain i perpendicular to the x and y axis. Futhermore, the 
resulting force supported by the domain i on the vertical axis z is equal to vYtl 2 . Therefore, the resulting 
elastic energy related to domain i and induced by the lattice mismatch e is given by the work of this elastic 
force, for a characteristic displacement el in the horizontal plane : 

E = Ye(2£hi + vl 2 )t£ (3) 

Consequently, the energy of a domain i with regard to its neighbors j might be written : 

NN 

E i = £ 2 {^-2Ye 2 )Y / (hi-h j ) (4) 
i=i 



3 Discussion 

The formation of QD's without nucleation is achieved through the evolution of the roughness of a 
thin film until nanometer scale islands are clearly identified. In the numerical models, the evolution of 
the roughness of a film is simulated by the evolution of the heights of the discrete mesoscopic domains 
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composing the virtual thin film. The energy of a mesoscopic domain i of a thin film synthesized via sol-gel 
route is calculated with equation (1). One can see that the roughness of the surface of the film is driven 
by surface diffusion and surface tension. The energy of a mesoscopic domain i of a thin film synthesized 
via a deposition process is calculated through equation (4). This equation allows one to understand that the 
roughness of the film is driven by a competition effect between surface and stress energies. 

The probability of changing the height of a domain with regard to equation (1) follows a monotonic 
tendency with regard to the intrinsic parameters of the thin film. Figure (1) presents the evolution of the is- 
landing of numerical thin films deposited on a substrate by sol-gel dip-coating after 10 7 MCS as a function 
of the initial thickness of the film. The annealing temperature is fixed, so all the parameters of the equa- 
tion are fixed. One can see that as the initial thickness of the film increases, the substrate is less dewetted. 
Indeed, as the initial thickness increases (h — 4nm), the amplitude of the roughness which appears at the 
top of the layer can not reach the substrate as fast as for thiner films (h = lnm). 

The probability of changing the height of a domain with regard to equation (4) depends on the relative 
values of the parameters j s /h and 2Ye 2 . Indeed, if j s /h > 2Ye 2 then the energy of the domain i is 
expressed as Ei — K J2j=i (hi ~ hj) (where K — j s /h — 2Ye 2 ) with K > 0, whereas K < otherwise. 
On the one hand, if the lattice mismatch e is low and/or with a low film thickness (K > 0), then the effect of 
surface energy is predominant. Therefore the probability of roughening is low and the film remains almost 
flat during the growing of the film which is described through the Frank-van-der-Merke (F-M) deposition 
process (or ALD Atomic Layer Deposition process). On the other hand, if the numerical values of e and/or 
film thickness are high (K < 0), then the probability of roughening increases ; therefore the effect of stress 
field inside the film is predominant which leads to the formation of islands. The surface roughening stage 
can either initiate from the top of a wetting layer, which is called the Stranski-Krastanov (S-K) growth 
mode, or directly from the top of the substrate without previous deposition of a wetting layer (e very high) 
which is the Volmer- Weber (V-W) growth mode. Figure (2) illustrates the numerical calculations of the 
islanding process through the effect of the initial film thickness (top, e=0.04) while the bottom presents 
this evolution with regard to lattice mismatch (bottom, h=2nm) after 10 7 MCS. We focus on the transition 
of K from positive to negative values as the islands are formed when K < 0. One can note that the increase 
of h or e leads to a higher number of QD's, which is contrary to the evolution of a thin film synthesized via 
sol-gel processing where the higher thickness of the film leads to a lower number of QD's. Moreover, it is 
clear that QD's formation is more sensitive to the variation of the lattice mismatch than to the variation of 
the initial thickness of the film as Ei is a function of e 2 and 1/h. Thus the model describes qualitatively 
the evolution from a F-M or S-K growth mode (K > and K w respectively) to a V-W growth mode 
K < as the number of QD's synthesized grows with the increase of e and/or h. Those tendencies have 
already been reported for the Ge/Si system [ 18 1. 

4 Conclusion 

Our simple energetic models allow us to describe the main tendency of QD's formation from thin film 
synthesized either via deposition processes (PVD-CVD) or sol-gel dip-coating. 

The main difference between those approaches is that thin films synthesized via sol-gel route lead 
to a greater number of QD's when they are thin while thin films synthesized via deposition processes 
demonstrate the opposite behaviour. This result is explained through equations (1) and (4) which point 
out that surface roughening of thin films synthesized via sol-gel route is driven by surface diffusion and 
surface tension, whereas it is the consequence of a competitive mechanism between elastic stress energy 
and surface tension for thin films synthesized layer by layer in PVD or CVD processes. 
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Fig. 1 Simulation of the islanding process of numerical thin films without deposition (sol-gel) as a function of the 
initial thickness of the film after 10 7 MCS. 
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Fig. 2 Simulation of the islanding process of numerical thin films with deposition (PVD-CVD) as a function of the 
initial thickness of the film (top) and of the lattice mismatch (bottom) after 10 7 MCS. 
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Table 2: Upright greek letters 
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Table 3: Boldface variants of slanted greek letters 
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Table 4: Boldface variants of upright greek letters 
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